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We address dilepton production in Au+Au collisions at -^/sjvjv = 200 GeV by employing the 
parton-hadron-string dynamics (PHSD) off-shell transport approach. Within the PHSD one goes 
beyond the quasiparticle approximation by solving generalized transport equations on the basis of 
the off-shell Kadanoff-Baym equations for the Green's functions in the phase-space representation. 
The approach consistently describes the full evolution of a relativistic heavy-ion collision from the 
initial hard scatterings and string formation through the dynamical deconfinement phase transition 
to the quark-gluon plasma (QGP) as well as hadronization and to the subsequent interactions in 
the hadronic phase. With partons described in the PHSD by the dynamical quasiparticle model 
(DQPM) - matched to reproduce lattice QCD results in thermodynamic equilibrium - we calculate, 
in particular, the dilepton radiation from partonic interactions through the reactions qq — > 7*, 
99 7* + 9 an d 19 7*9 (<7<? ~~ ► 7*?) m the early stage of relativistic heavy-ion collisions. By 
comparing our results to the data from the PHENIX Collaboration, we study the relative importance 
of different dilepton production mechanisms and point out the regions in phase space where partonic 
channels are dominant. Furthermore, explicit predictions are presented for dileptons within the 
acceptance of the STAR detector system and compared to the preliminary data. 

PACS numbers: 25.75.-q, 25.75.Cj, 25.75.Nq, 24.85,+p 



I. INTRODUCTION 

The nature of confinement and the phase transition 
from a partonic system of quarks, antiquarks and gluons, 
the so-called quark-gluon plasma (QGP), to interacting 
hadrons, as occurring in relativistic nucleus-nucleus col- 
lisions, is a central topic of modern high energy physics. 
Already some decades ago dileptons were suggested as 
useful probes of the properties of the QGP [iHal- Since 
dileptons arc emitted over the entire history of the heavy- 
ion collision, from the initial nuclcon-nucleon collisions 
through the hot and dense (partonic) phase and to the 
hadron decays after freeze-out, microscopic transport 
models have to be applied for disentangling the various 
sources that contribute to the final dilepton spectra seen 
in experiments. 

In the present work, we study dilepton production 
in Au+Au collisions at ^/JWn = 200 GeV within the 
PHSD off-shell transport approach by including the col- 
lisional broadening of vector mesons, microscopic sec- 
ondary multi-meson channels, and the radiation from the 
strongly interacting QGP through the interactions of dy- 
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namical quasiparticles having broad spectral functions in 
line with the degrees of freedom propagated in the trans- 
port approach. 

The PHENIX Collaboration has presented dilep- 
ton data from pp and Au+Au collisions at the 
Relativistic-Heavy-Ion-Collidcr (RHIC) energy of 
y / s]vlv=200 GeV [B-dl , which show a large enhancement 
relative to the scaled pp collisions in the invariant mass 
regime from 0.15 to 0.6 GeV/c 2 . The question that we 
aim at investigating in this work is whether this excess 
can be attributed to the gluon Compton scattering or 
other partonic interaction processes, which due to the 
two-particle finite state contribute at low masses as well. 

Moreover, if a realistic partial loss of the D- and 
-D-meson correlations due to their rescattering is taken 
into account, one has to conclude from the data of the 
PHENIX Collaboration that there exists another domain 
of invariant masses, in which the measured dilepton yield 
in Au+Au collisions is underestimated by the scaled 
yield from the p + p collisions, i.e. at masses from 1 
to4GeV/c 2 [1G3. The microscopic calculations within 
the parton-hadron transport approach here will answer 
if this discrepancy can be accounted for by the partonic 
sources of dileptons. 

In Ref. [Tl[ we have recently presented results from the 
PHSD for the dilepton spectrum produced in In+In col- 
lisions at 158 A-GeV and compared them to the NA60 
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data [l2j, ll3( . We recall that by employing the hadron- 
string dynamics (HSD) transport approach that does not 
include explicit partonic contributions to low mass dilep- 
ton production in relativistic heavy-ion collisions, it was 
shown in Ref. [l4| that the NA60 data for the invari- 
ant mass spectra of pairs from In+In collisions at 
158 A-GcV indicated an in-medium modification of the p- 
meson according to the 'melting' scenario [12|. The more 
recent PHSD calculations in Ref. [Tl[ , which took into ac- 
count the phase transition to and the radiation from the 
partonic phase, confirmed the HSD results that the spec- 
trum at invariant masses in the vicinity of the p peak was 
better reproduced by the p meson yield if a broadening 
of its spectral function in the medium was assumed. On 
the other hand, the spectrum at M > 1 GeV was shown 
to be dominated by partonic sources. Moreover, the in- 
clusion of the partonic dilepton sources made it possible 
to reproduce in the PHSD the effective temperature or 
the inverse slope parameter of the transverse momentum 
spectrum of dileptons in the intermediate mass region. 
Furthermore, for dileptons of low masses (M < 0.6 GeV), 
a sizable contribution from partonic processes, particu- 
larly the quark annihilation with gluon brcmsstrahlung 
in the final state, was found, and this provides another 
possible window for probing the properties of the sQGP. 

In the present work, our previous findings will be re- 
examined by comparison to the dilepton measurements 
at RHIC energies. Previous HSD results for e + e~ pairs 
in Au+Au collisions in comparison to the data from the 
PHENIX Collaboration @,||] were presented in Ref. [l4j . 
Whereas the total yield is quite well described in the low- 
mass region from the pion Dalitz decay as well as around 
the u, <fi an d J/ty masses, the HSD clearly underesti- 
mates the measured spectra in the regime from 0.2 to 0.6 
GeV by approximately a factor of 5 for central Au+Au 
collisions. After including the in-medium modification 
of vector mesons, we obtain a total spectrum which is 
only slightly enhanced compared to the 'free' scenario 
of using the vector meson properties in the vacuum (see 



Fig. 6 of Ref. |1J]). The low mass dilepton spectra from 



Au+Au collisions obtained by the PHENIX Collabora- 
tion at RHIC are thus clearly underestimated in the in- 
variant mass region of 0.2 to 0.6 GeV in the 'collisional 
broadening' scenario as well as in the 'dropping mass + 
collisional broadening' model. We mention that the HSD 
results for the hadronic production of low mass dileptons 
are very close to those calculated by van Hees and Rapp 
as well as by Dusling and Zahed [15[ (cf. the comparison 
in Refs. HG3). 

At higher masses from 1 to 4 GeV, the dominant 
hadronic sources of lepton pairs are from the semi- 
leptonic decays of correlated D-mesons and the dilepton 
decays of charmonia. In estimating this contribution, the 
experimental information was used in Ref. [Io| to deter- 
mine the yields of charmed hadrons. Additionally, the 
effect of D meson rescattering in the hadronic matter 
was included in the HSD transport approach [l?], EH to 
estimate the probability of surviving correlated semilcp- 



tonic charm decays in heavy-ion collisions. The resulting 
dilepton yield between the cj> and J/4" peaks - after in- 
cluding the semi-leptonic decays of correlated D-mesons 
in the HSD - was found to underestimate the PHENIX 
data by approximately a factor of two [Toj . 



Another open question to be answered by the micro- 
scopic transport calculations is the determination of 'win- 
dows' in phase space for observing dileptons from the 
quark-gluon plasma that possibly overshine the hadronic 
sources. It has been originally suggested that a substan- 
tial thermal yield from the deconfined phase might be 
seen in the invariant mass region between the <fi and J7\E' 
peaks Q, while the spectrum at lower masses was dom- 
inated by meson decays. On the other hand, the calcu- 
lations in Refs. [ll|, [lj| Ho| pointed to a possible second 
region of phase space for the observation of the ther- 
mal QGP source at masses ~ 0.3 — 0.6 GeV. In order to 
clarify whether these dileptons of masses 0.3 — 0.6 GeV 
can be observed among the background of dileptons from 
hadronic decays, a study within a transport approach 
that incorporates dilepton production from the (non- 
equilibrium) partonic phase, hadronic decays and the mi- 
croscopic secondary hadronic interactions - including the 
"47r" channels - thus appears appropriate. 



The PHSD [2l|, |22( transport approach, which incor- 
porates the relevant off-shell dynamics of vector mesons 
and the explicit partonic phase in the early hot and dense 
reaction region as well as the dynamics of hadronization, 
allows for a microscopic study of various dilepton pro- 
duction channels in non-equilibrium matter. The PHSD 
off-shell transport approach is particularly suitable for in- 
vestigating the enhanced production of lepton pairs in the 
invariant mass range 0.3 < M < 0.7 GeV/c 2 that is seen 
in experiments, since it incorporates various scenarios for 
the modification of vector mesons in a hot and dense 
medium. In the present work, we calculate the spectra of 
dileptons produced in the course of Au+Au collisions at 
\/snn— 200 GeV from the partonic and hadronic sources 
by including the partonic channels and the multi-meson 
channels besides the usual hadron decay channels. By 
consistently treating in the same microscopic transport 
framework both partonic and hadronic phases of the col- 
liding system, we are aiming to determine the relative 
importance of different dilepton production mechanisms 
and to point out the regions in phase space where par- 
tonic channels are dominant. 



The paper is organized as follows. In Sec. [Til we gi ye 
a brief description of the PHSD approach as well as the 
hadronic and partonic sources of dilepton production in- 
corporated in the PHSD. We then compare in Sec. lIIII rc- 
sults from our calculations to the available experimental 
data for Au+Au collisions at top RHIC energy. Finally, 
our conclusions are presented in Sec. IIV1 
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II. THE PHSD TRANSPORT APPROACH 

The PHSD (HHH is an off-shell transport model that 
consistently describes the full evolution of a rclativistic 
heavy-ion collision from the initial hard scatterings and 
string formation through the dynamical deconfinement 
phase transition to the quark-gluon plasma as well as 
hadronization and to the subsequent interactions in the 
hadronic phase. In the hadronic sector, the PHSD is 
equivalent to the HSD transport approach [23Tj25j that 
has been used for the description of pA and AA collisions 
from SIS to RHIC energies and has lead to a fair repro- 
duction of measured hadron abundances, rapidity distri- 
butions and transverse momentum spectra. In particular, 
as in the HSD, the PHSD incorporates off-shell dynamics 
for vector mesons [26[ and a set of vector-meson spectral 
functions [2lJ that covers possible scenarios for their in- 
medium modifications. The transition from the partonic 
to hadronic degrees of freedom is described by covariant 
transition rates for the fusion of quark-antiquark pairs 
to mesonic resonances or three quarks (antiquarks) to 
baryonic states, i.e., the dynamical hadronization j28| . 
Note that due to the off-shell nature of both partons and 
hadrons, the hadronization process obeys all conservation 
laws (i.e., the 4-momentum conservation and the flavor 
current conservation) in each event, the detailed balance 
relations, and the increase in the total entropy S. The 
transport theoretical description of quarks and gluons in 
the PHSD is based on a dynamical quasiparticlc model 
(DQPM) for partons that is matched to reproduce the 
lattice QCD (1QCD) results for a quark-gluon plasma in 
thermodynamic equilibrium. The DQPM provides the 
mean-fields for gluons/quarks and their effective 2-body 
interactions in the PHSD. For details about the DQPM 
model and the off-shell transport we refer the reader to 
Ref. [29[. We stress that a non- vanishing width 7 in the 
partonic spectral functions is the main difference between 
the DQPM and conventional quasiparticle models [30| . 
Its influence on the collision dynamics is essentially seen 
in the correlation functions, e.g., in the stationary limit, 
the correlation involving the off-diagonal elements of the 
energy-momentum tensor T kl defines the shear viscosity 
?y of the medium 31| . Here a sizeable width is mandatory 
to obtain a small ratio of the shear viscosity to entropy 
density 77/s, which results in a roughly hydrodynamical 
evolution of the partonic system in PHSD (28|. The fi- 
nite width leads to two-particle correlations, which are 
taken into account by means of the generalized off-shell 
transport equations [26| , that go beyond the mean field 
or Boltzmann approximation [29l . |32| . 



A. Partonic sources of dileptons in PHSD 

In the scope of the one- and two-particle interactions, 
dilcpton radiation by the constituents of the strongly in- 
teracting QGP proceeds via following elementary pro- 
cesses: the basic Born q + q annihilation mechanism, 



gluon Compton scattering (q + g —> 7* + q and q + g —> 
7* + q), and quark and anti-quark annihilation with the 
gluon Bremsstrahlung in the final state (q + q — > g + 7*). 
In the on-shcll approximation, one would use perturba- 
tive QCD cross sections for the processes listed above. 
However, in the strongly interacting QGP the gluon 
and quark propagators differ significantly from the non- 
interacting propagators. Accordingly, we have calculated 
in Rcfs. [33, [34| the off-shell cross sections for dilcpton 
production in the partonic channels by off-shell partons, 
using the phenomenological parametrizations (from the 
DQPM) for the quark and gluon propagators and their 
interaction strength. 

In Refs. [Hldll it was shown that the finite quark and 
gluon masses modify the magnitude as well as the M— 
and pt~ dependence of the cross sections of the above 
mentioned processes compared to the perturbative re- 
sults for massless partons (cf. Figs. 3 and 4 of Ref. [H|). 
The modifications are large at lower M 2 and at the edges 
of the phase space. It was shown that the most promi- 
nent effect of the quark masses on the dimuon produc- 
tion cross sections in the Born mechanism (q + q — > 7*) 
was a sharp threshold value for the invariant mass of the 
dilepton pair M m i n = mi + 7712. On the other hand, 
the finite masses of the quark and antiquark produce 
additional higher-twist corrections to the cross section, 
which decrease with increasing M 2 , so that the off-shell 
cross sections approach the leading twist - on-shell - re- 
sult in the limit of high dilepton masses. In Fig. 4 of 
Ref. [34| , an analogous comparison for the 2 — > 2 process 
q + q — > 7* + g was shown by plotting the off-shell (i.e. 
with finite masses for the quarks and gluons) cross section 
for the quark annihilation with gluon bremsstrahlung in 
the final state of various values of the quark and gluon 
off-shellncss (masses) and the corresponding on-shell re- 
sult. As found in Ref. 34 1. the maximum pair mass shifts 
to a lower value as a result of producing a massive gluon 
in the final state. For the rest of the M values, the ef- 
fect of the quark and gluon masses is about 50%. For 
m q /g — > 0, the cross section approaches the leading twist 
pQCD result. 

The question of the effect of a finite parton width - 
which parametrizes their interaction rate and correlation 
as well as multiple scattering - on dilepton rates in heavy- 
ion collisions has been addressed in Ref. [HI by convo- 
luting the off-shell cross sections with phenomenological 
spectral functions A(m q ) and A(rn g ) for the quarks and 
gluons in the quark-gluon plasma, respectively, and with 
parton distributions in a heavy-ion collision. The finite 
width of the quasiparticles was found to have a sizable 
effect on the dilepton production rates. In particular, 
the threshold of the Drell-Yan contribution was "washed 
out" . Also, the shape and magnitude of the 2 — > 2 pro- 
cesses (q + q — > <7 + 7* and q + g 3 + 7*) were modified. 
One further observed that the contribution of the gluon 
Compton process q + g — > q + 7* to the rates was small 
compared to that of q + q annihilations. 

In Ref. [ll[ , we have implemented the cross sections ob- 



4 



taincd in Rcfs. |33j,l34| into the PHSD transport approach 
in the following way: Whenever the quark-antiquark, 
quark-gluon and antiquark-gluon collisions occur in the 
course of the Monte-Carlo simulation of the partonic 
phase in the PHSD, a dilepton pair can be produced ac- 
cording to the off-shell cross sections which depend, 
in addition to the virtualitics of the partons involved, on 
the energy density in the local cell where the collision 
takes place. The local energy density governs the widths 
of the quark and gluon spectral functions as well as the 
strong coupling in line with the DQPM. 



B. Hadronic sources of dileptons in PHSD 

In the hadronic sector, the PHSD is equivalent to the 
HSD transport approach [23l - l25| . The implementation of 
the hadronic decays into dileptons (it-, T}-, rj'-, uj-, A-, ai- 
Dalitz, p -> l + l~, uj -> l + l~ , and cf> -> l + l~) in HSD (and 
PHSD) is described in detail in Refs. 0, S3- For the 
treatment of the leptonic decays of open charm mesons 
and charmonia, we refer the reader to Refs. [HI Ell 13 
for a detailed description. 

The PHSD off-shell transport approach is particularly 
suitable for investigating the different scenarios for the 
modification of vector mesons in a hot and dense medium. 
Just as the HSD model, the PHSD approach incorpo- 
rates the off-shell propagation for vector mesons as de- 
scribed in Ref. [26]. In the off-shell transport, the hadron 
spectral functions change dynamically during the propa- 
gation through the medium and evolve towards the on- 
shell spectral functions in the vacuum. As demonstrated 
in Ref. [13], the off-shell dynamics is important for res- 
onances with a rather long lifetime in the vacuum but 
strongly decreasing lifetime in the nuclear medium (es- 
pecially uj and <p mesons) and also proves vital for the 
correct description of dilepton decays of p mesons with 
masses close to the two pion decay threshold. For a de- 
tailed description of the off-shell dynamics and the im- 
plementation of vector-meson modifications in medium, 
we refer the reader to Refs. [H GJ, SI, S3, S3 • 

In Ref. [ll|, the extension of the hadronic sources in 
the PHSD for dilepton production to include secondary 
multi-meson interactions by incorporating the channels 
ttuj — > l + l~ , nai — > l + l~ , and pp — > l + l~ is described in 
detail. These so-called '47r channels' for dilepton produc- 
tion are incorporated in the PHSD on a microscopic level 
rather than assuming thermal dilepton production rates 
and incorporating a parametrization for the inverse reac- 
tion p + + p~ — > Air' s by employing the detailed balance 
as in Refs. [H, S3- By studying the electromagnetic 
emissivity (in the dilepton channel) of the hot hadron 
gas, it was shown in Refs. [13, SH that the dominating 
hadronic reactions contributing to the dilepton yield at 
the invariant masses above the <f> peak are the two-body 
reactions of 7r + p, tt + uj, p + p, and tt + ai . This conclu- 
sion was supported by the subsequent study in a hadronic 
relativistic transport model [39j . Therefore, we have im- 



plemented the above listed two-meson dilepton produc- 
tion channels in the PHSD approach. In addition, some 
higher vector mesons (p' etc.) are tacitly included by us- 
ing phenomenological form factors that are adjusted to 
the experimental data (cf. Ref. [ll| ) ■ Specifically, we de- 
termined the cross sections for the mesonic interactions 
with dileptons in the final state using an effective La- 
grangian approach, following the works of Refs. [37l[39j. 
In order to fix the form factors in the cross sections for 
dilepton production by the interactions of tt + p, tt + to, 
p + p and Trai, we used the measurements in the detailed- 
balance related channels: e + e~ — > tt + p, e + e~ — > tt + uj, 
e + e~ — > p+p, and e + e~ — > TT+a\. Note that we fitted the 
form factors while taking into account the widths of the p 
and a\ mesons in the final state by convoluting the cross 
sections with the (vacuum) spectral functions of these 
mesons in line with Ref. [40 ] (using the paramctrizations 
of the spectral functions as implemented in the HSD and 
described in [4l|). In Fig. 5 of Ref. [ll| we have presented 
the resulting cross sections, which are implemented in the 
PHSD. 



III. COMPARISON TO DATA 

Let us first note that the bulk properties of heavy- 
ion collisions at the top RHIC energy, such as the num- 
ber of charged particles as well as their rapidity, px, v 2 
and transverse energy distributions, are rather well de- 
scribed by the PHSD; and we refer to Ref. [22j for an 
extended and detailed comparison to the data from the 
PHOBOS, STAR and PHENIX collaborations. Since 
the 1QCD equation of state employed in the PHSD has 
a crossover transition, the PHSD calculations show a 
rather long QGP phase in central Au+Au collisions at 
y/sl^N = 200 GeV (cf. Fig. 6 of Ref. [13]) with the par- 



tonic degrees of freedom dominating for about 5-7 fm/c. 
We recall that dilepton production in the elementary pp 
channel is also well under control in the PHSD as was 
previously demonstrated in Ref. [l4[ . 

In Fig. fl] we show our results for the invariant mass 
spectra of inclusive dileptons in Au+Au collisions in the 
invariant mass region M = — 1.2 GeV for the accep- 
tance cuts on single electron transverse momenta p e T, 
pseudorapiditics rj e , azimutal angle <j) e , and dilepton pair 
rapidity y, 

PeT > 0.2 GeV, 
|r? e | < 0.35, 

-3tt/16 < (j) e < 5tt/16, IItt/16 < <j) e < 19tt/16, 

\y\ < 0.35. (1) 

In this region, the dilepton yield in the PHSD is domi- 
nated by hadronic sources and essentially coincides with 
the earlier HSD result [l4[. There is a sizeable discrep- 
ancy between the PHSD calculations and the data from 
the PHENIX Collaboration in the region of masses from 
0.2 to 0.6 GeV. The discrepancy is not amended by the 
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FIG. 1: The PHSD results for the invariant mass spectra 
of inclusive dileptons in Au+Au collisions at ^/sjvjv = 200 
GeV within the PHENIX acceptance cuts [Eq.Q] as given in 
the main text in comparison to the data from the PHENIX 
Collaboration 0, @] for invariant masses M = — 1.2 GeV. 
The different lines indicate the contributions from different 
channels as specified in the figure. 



inclusion of the radiation from the QGP as well as from 
correlated D-meson decays, since the latter contributions 
are 'over-shone' by the radiation from hadrons integrated 
over the entire evolution of the collision. 

In contrast, the partonic radiation as well as the yield 
from correlated D-meson decays are dominant in the 
mass region M = 1 — 4 GeV as seen in Fig. [21 i.e., in 
the mass region between the <fi and J/ty peaks. The 
dileptons generated by the quark- ant iquark annihilation 
in the sQGP constitute about half of the observed yield 
in this intermediate-mass range. For M > 2.5 GeV the 
partonic yield even dominates over the D-meson contri- 
bution. Thus, the inclusion of the partonic radiation in 
the PHSD fills up the gap between the hadronic model 
results [H [HI and the data of the PHENIX Collabora- 
tion for M > 1 GeV. Note that the collisional broadening 
scenario for the modification of the p-meson was used in 
the calculations presented in Figs. [T] and [5] 



A. Centrality dependence 

Next we investigate the centrality dependence of dilep- 
ton production in heavy-ion collisions at the top RHIC 
energy for the centrality cuts specified in Eq.([T|). While 
results from the PHSD calculations are in a reasonable 
agreement with the PHENIX data at 20-40% centrality, 
the data show an increasing excess over the PHSD re- 
sults for the 10-20% centrality in the mass regime from 
0.15-0.7 GeV, and this excess becomes even more dra- 
matic for the most central (0-10% centrality) collisions 
(cf. Figs. 3(a)|3(c)" ). For all these centrality cuts, the 
contribution of partonic channels from the PHSD is sub- 
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FIG. 2: Same as Fig. [T] for invariant masses M = — 4 GeV. 



leading for 0.15 GeV < M < 0.7 GeV and cannot be 
responsible for the excess dileptons seen experimentally 
by the PHENIX Collaboration. In short, the early ex- 
pectation of a partonic signal in the low mass dilepton 
spectrum is not verified by the microscopic PHSD calcu- 
lations. 



B. Channel decomposition in central collisions 



In order to elucidate the relative importance of the dif- 
ferent hadronic sources of the excess dileptons in heavy- 
ion collisions at the top RHIC energy, we show in Fig. 0] 
the channel decomposition of the main hadronic contri- 
butions to the dilepton rates in central Au+Au collisions 
at y/SNN = 200 GeV without any cuts on dilepton mo- 
menta and rapidity and without including the finite mass 
resolution of the PHENIX detector system. Here the 
blue line with its sharp peaks from the decay of vector 
mesons and a smooth background from the Dalitz de- 
cays of hadrons becomes rather insignificant above the 
</>-meson mass. Also the 47r-channels (tt + a\ , n + u>, -k + p, 
and p + p) are clearly subleading in the intermediate- 
mass region. Here the partonic channels - dominated by 
q + q ^ e + e~ - constitute about half of the dilepton yield 
and have about the same contribution as that from cor- 
related D-meson decays (not shown in Fig. - Note that 
the contribution from partonic channels is approximately 
exponential in mass for 1 GeV < M < 2.5 GeV and might 
be interpreted as being due to 'thermal radiation' from 
the sQGP. However, the PHSD calculations do not indi- 
cate that a thermal equilibrium has been achieved on the 
partonic level. 
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FIG. 3: The PHSD results for the dilepton invariant mass spectra in Au+Au collisions of different centralities at ^/snn = 200 
GeV in comparison to the data from the PHENIX Collaboration 0, H| . 



C. Transeverse momentum distributions 

The PHENIX Collaboration has also accessed the 
information on the transverse momentum dependence 
of dilepton production by measuring the py-spectra of 
dilcptons in different bins of invariant mass M. In 
Fig. [5]wc show the measured transverse momentum spec- 
tra of dilcptons for minimum bias Au+Au collisions at 
•J snn = 200 GcV (symbols) in comparison with the 
spectra from the PHSD (lines) for six mass bins as indi- 
cated in the figure. Whereas the PHSD can well describe 
the dilepton spectra in the mass intervals [0,100 McV] 
and [810 MeV, 990 McV], it underestimates the low p T 
dileptons in the other mass bins, particularly in the mass 
bins [300 MeV, 500 MeV]. On the other hand, high p T 
dilcptons are reproduced quite well by the PHSD calcu- 
lations. Wc conclude that the missing dilepton yield for 
masses from 0.15 to 0.6 GeV is essentially due to a se- 
vere underestimation of the data at low px by up to an 
order of magnitude. We recall that at top SPS energies 
the low pt dilepton yield could be attributed to tt — tt 
annihilation channels, i.e. to the soft hadronic reactions 
in the expansion phase of the system. These channels 



are, however, insufficient to describe the very low slope 
of the pt spectra at the top RHIC energy 



D. Comparison with other models 

Since the authors have worked on the topic of dilepton 
production within various approaches, it is instructive to 
discuss the results especially in the low mass regime for 
the same centralities and within the acceptance of the 
PHENIX detector system. In Fig. [5] we thus present the 
invariant mass spectra of inclusive dileptons in Au+Au 
collisions at i/sjvjv = 200 GeV in the low mass region 
(M = 0-1.2 GeV) as calculated within 1) PHSD (red solid 
line) taking into account the in-mcdium modification of 
the p as well as the dilepton radiation from the partonic 
phase; 2) HSD in the free-p scenario; and 3) the extended 
thermal model [loj in comparison to the data from the 
PHENIX Collaboration @,t|. 

The extended statistical hadronization model 
(ESHM) Q3 is an extension to the statistical hadroniza- 
tion model (SHM) which has been applied [42T(55j 
to high-energy elementary and especially heavy-ion 
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FIG. 4: Dilepton invariant mass spectra from central Au+Au 
collisions at ^Jsnn = 200 GeV integrated over rapidity and 
transverse momenta as calculated in PHSD. The contribu- 
tions from multi-pion channels, quark-antiquark annihilation 
and hadron decays into dileptons are shown separately. 
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FIG. 5: The PHSD results for the transverse-momentum 
spectra of dileptons from minimum bias Au+Aucollisions at 
v / ijviv = 200 GeV in different mass bins compared to the data 
from the PHENIX Collaboration 



collision experiments in order to calculate the yields 
of different hadron species. In the SHM, the state of 
the "thermal" fireball is specified by its temperature T, 
volume V and the chemical potentials [1b, (J-q, and /is 
for baryon, electric and strangeness charges, respectively. 
While /is and (iq are zero in central Au+Au collisions 
at y/sjjjf = 200 GeV, /i B is about 30 MeV. The effect of 
the strangeness under-saturation parameter or fugacity 
7s on the di-clectron invariant mass spectrum as a 
function of centrality was studied in detail in Ref. [Tfj| 
and the effect was found to be moderate. We employ 
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FIG. 6: Invariant mass spectra of inclusive dileptons in 
Au+Au collisions at y / ijvjv = 200 GeV in the low mass region 
(M=0-1.2 GeV) as calculated in three models - PHSD tak- 
ing into account in-medium modification of rho mesons and 
dilepton radiation from the partonic phase (solid line); HSD 
in the free-rho scenario (dashed line) ; and the extended ther- 
mal model [ro(l - compared to the data from the PHENIX 
Collaboration^, ■ 



the value 7s=0.6 in this work for the min bias Au+Au 
collisions. The overall normalization (fireball volume) at 
different centralities was fitted to experimental data in 
Ref. [l(| and we use the same values throughout. For 
the temperature we use the value T = 170 MeV. 

Since the measured rapidity and transverse momentum 
spectra of hadrons emitted in the high energy collision 
experiments do not resemble thermal distributions, the 
SHM has been extended in Ref. by boosting (event by 
event) the "fireball" along the beam axis so that the ra- 
pidity distributions of pions become compatible with the 
BRAHMS measurements (5(|. Also, the problem that 
the SHM tends to over-populate the low pr part of the 
spectrum compared with the experimental distributions 
was solved in the ESHM by assuming that the created 
clusters' transverse momentum is normally distributed 
with the width fitted together with the system volume V 
to the PHENIX data [531 m -P + P collisions and in 11 dif- 
ferent centrality classes in the case of Au+Au collisions. 
For further details we refer the reader to Ref. [loj . 

We find in Fig.[B]that the HSD and the ESHM give ap- 
proximately the same results on the level of 30% for the 
dilepton invariant mass spectra. This might be surprising 
since the HSD includes not only the direct and Dalitz de- 
cays of hadrons but also meson-meson and meson-baryon 
channels for dileptons. Indeed, the enhancement of the 
HSD result from 0.55 to 0.75 GeV can be traced back to 
pion-pion annihilation which, however, gives only a small 
contribution at the top RHIC energy. Our actual PHSD 
calculations show some more dilepton yield in the p-mass 
regime as a result of the broadened p spectral function 
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FIG. 7: The PHSD results for the invariant mass spectra of 
inclusive dileptons in Au+Au collisions at ^snn = 200 GeV 
for M = 0-1.2 GeV and - 80 % centrality within the cuts 
of the STAR experiment, see [Eq.(|2j] in the main text. The 
preliminar y d ata from the STAR Collaboration are adopted 
from Ref. 16011. 




M [GeV/c ] 

FIG. 8: Same as Fig. [7] for M=0-4 GeV. 



employed in the calculations. In the 'free /?' scenario, the 
results from HSD and PHSD are identical within statis- 
tics since the partonic channels give only a small con- 
tribution in this mass range. The conclusion that the 
dilepton spectrum at masses below 1 GeV is dominated 
by the hadronic sources is also supported by the studies 
in other available models [H, Hil • 



E. Predictions for STAR and comparison to 
preliminary data 

The PHSD calculations allow us to match with the 
different experimental conditions and thus to provide a 
'theoretical link' between the different measurements. To 



.4 0.6 2 0.8 1.0 1 

M [GeV/c ] 

FIG. 9: The PHSD results for the invariant mass spectra of 
central dileptons in Au+Au collisions at ^/snn = 200 GeV 
for M = 0-1.2 GeV and - 10 % centrality within the cuts 
of the STAR experiment, see Eq.[((2])] in the main text. The 
preliminar y d ata from the STAR Collaboration are adopted 
from Ref. [6(| . 




FIG. 10: Same as Fig. [9] for M = - 4 GeV. 



this extent, we have provided the differential data-tables 
for our theoretical predictions on the web site [6l[ so that 
any acceptance cuts and experimental mass and trans- 
verse momentum resolutions can be applied. 

The STAR Collaboration at RHIC has recently mea- 
sured dileptons from Au+Au collisions at tJsnn = 200 
GeV with the acceptance cuts on single electron trans- 
verse momenta p e r, single electron pseudorapidities r] e 
and the dilepton pair rapidity y 

0.2 < p eT < 5 GeV, 

\Ve\ < 1, 

\y\ < i- (2) 

Our predictions for the dilepton yield within these cuts 
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are shown in Figs. [7] and [8] for 0-80% ccntrality and in 
Figs . 151 and [TU1 for 0-10% centrality. One can observe gen- 
erally a good agreement with the preliminary data from 
the STAR Collaboration [goj] for - 80% centrality in the 
whole mass regime. Surprisingly our calculations are also 
roughly in line with the low mass dilepton spectrum from 
STAR in case of central collisions whereas the PHSD re- 
sults severely underestimate the PHENIX data for the 
cuts given in Eq.(l) (cf. Fig. 3). The observed yield 
from STAR can be accounted for by the known hadronic 
sources, i.e. the decays of the ir , rj, rf , u, p, (p and a\ 
mesons, of the A particle and the semileptonic decays of 
the D and D mesons, where the collisional broadening of 
the p meson is taken into account! At first sight this ob- 
servation might point towards an inconsistency between 
the data sets from PHENIX and STAR but we have to 
stress that the actual experimental acceptance cuts are 
more sophisticated than those given in Eqs.(l) and (2). 
This problem will have to be investigated closer by the 
experimental collaborations. Furthermore, the upgrade 
of the PHENIX experiment with a hadron blind detec- 
tor [62| should provide decisive information on the origin 
of the low mass dilcptons produced in the heavy-ion col- 
lisions at yfs = 200 GcV. 

We also observe a slight overestimation of the dilepton 
yield from PHSD in 0-10% central collisions at masses 
from 1.3 to 1.8 GeV where the dominant contributions 
to the spectrum are the radiation from the sQGP and the 
semileptonic decays of the D and D mesons. We spec- 
ulate that the suppression of dilcptons from the D and 
D mesons might be underestimated in the PHSD calcu- 
lations in central collisions. The upgrade of the STAR 
detector [63| will be promising in independently measur- 
ing the correlated D and D meson contributions. 



IV. SUMMARY 

In this study, we have addressed dilepton production in 
Au+Au collisions at yfsNN = 200 GeV by employing the 
Parton-Hadron-String Dynamics (PHSD) off-shell trans- 
port approach. This work is a continuation of our ear- 
lier studies for heavy ion collisions at the SIS energies 
of 1-2 A GeVjp and the SPS energies from 40 to 158 
A GeV [ill, Il4l| essentially within the same dynamical 
transport model. Within the PHSD one solves gener- 
alized transport equations on the basis of the off-shell 
Kadanoff-Baym equations for effective Green's functions 
in phase-space representation (beyond the quasiparticlc 
approximation) for quarks, antiquarks and gluons as well 
as for the hadrons and their excited states. The PHSD 
approach consistently describes the full evolution of a 
relativistic heavy-ion collision from the initial hard scat- 
terings and string formation through the dynamical de- 
confinement phase transition to the quark-gluon plasma 
(QGP) as well as hadronization and the subsequent in- 
teractions in the hadronic phase. It was shown in previ- 
ous studies that the PHSD approach well describes the 



various hadron abundancies, their longitudinal rapidity 
distributions as well as transverse momentum distribu- 
tions from lower SPS to top RHIC energies [HI, [HI- Als o 
the collective flow ^(p*) is roughly in accordance with 
the experimental observations by the PHOBOS, STAR 
and PHENIX collaborations at RHIC 0. The latter 
findings allow to explore the dynamics of subleading/rarc 
probes within the dynamical environment of partons and 
hadrons during the complex time evolution of a relativis- 
tic heavy-ion collision. 

The present study has been devoted particularly to 
the calculation of dilepton radiation from partonic inter- 
actions through the reactions qq — > 7*, qq — > 7* + g and 
9.9 7*9 (93 — ► 7*9) m the early stage of relativistic 
heavy-ion collisions at the top RHIC energy. We recall 
that the differential cross sections for electromagnetic ra- 
diation have been calculated with the same propagators 
as those incorporated in the PHSD transport approach. 
By comparing our calculated results to the data from 
the PHENIX Collaboration, wc have studied the rela- 
tive importance of different dilepton production mecha- 
nisms and addressed in particular the 'PHENIX puzzle' 
of a large enhancement of dilcptons in the mass range 
from 0.15 to 0.6 GeV as compared to the emission of 
hadronic states. Our studies have demonstrated that the 
observed excess in the low mass dilepton regime cannot 
be attributed to partonic productions as expected ear- 
lier. Thus the 'PHENIX puzzle' still has no explanation 
from the theoretical approaches so far. The PHENIX 
enhancement is essentially due to dileptons of low trans- 
verse momentum in the mass range from 0.15 to 0.6 GeV 
and up to date finds no explanation by hadronic nor par- 
tonic reaction channels in PHSD that occur on top of 
the interactions in pp collisions during the evolution of 
relativistic heavy-ion collisions. 

we find 



Similar to our findings at SPS energies [11 



that the partonic dilepton production channels are visi- 
ble in the intermediate-mass region between the (f> and 
J/\I/ peaks. Their contribution is about the same as 
the correlated background from £>-meson decays. Sur- 
prisingly, this contribution appears to be exponential in 
mass from 1 to 2.5 GeV so that an interpretation of 'ther- 
mal radiation from the sQGP' might appear appropriate. 
However, the PHSD dynamics shows that no kinetic equi- 
librium is achieved on the partonic level in heavy-ion col- 
lisions at top RHIC energies, and such an interpretation 
has to be considered with care. 

In view of the fact that within different (statistical and 
dynamical) models we have not been able to find an ex- 
planation for the low mass 'PHENIX puzzle', the solu- 
tion has to be relegated to the experimental side. In fact, 
the STAR collaboration has taken independent dilepton 
data for centralitics different from the PHENIX measure- 
ments and also with different detector acceptances. Our 
PHSD calculations allow to match with the different ex- 
perimental conditions and thus to provide a 'theoretical 
link' between the different measurements. To this ex- 



tent, we have provided our predictions [61[ for the con 
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ditions of the STAR experiment, which happen to be in 
a rough agreement with the preliminary data from the 
STAR Collaboration [6(j. This finding opens up new 
problems that will have to be addressed from the exper- 
imental side. The upgrade of the PHENIX experiment 
with a hadron blind detector [62j and the upgrade of the 
STAR detector 63] for independently measuring the cor- 
related D and D meson contributions appear mandatory 
to shed some further light on the present 'puzzles'. 
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